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In this contribution the new event generation framework SHERPA will be presented. It aims at the full simulation 
of events at current and future high-energy experiments, in particular the LHC. Some results related to the 
production of jets at the Tevatron will be discussed. 



1. INTRODUCTION 

The observation and interpretation of multi- 
particle, multi-jet final states will be in the cen- 
tre of the physics programme at the LHC. They 
serve as signals or backgrounds for new physics; 
as an example consider the production and decay 
of heavy SUSY particles. This shift of focus to- 
wards higher multiplicities translates directly into 
new challenges for the simulation of such events 
and necessitates the construction of new simula- 
tion tools; for an overview on currently available 
tools cf. pp. 

The multi-purpose event generator SHERPA [2] is 
one of these new tools, which are currently un- 
der construction. Others, namely PYTHIA7 
H], PYTHIA8 and HERWIG++ are complete 
rewrites of the well-established codes PYTHIA [?| 
E] and HERWIG |9I1()| . extending or improving 
their physics content. To exemplify this, con- 
sider the new parton showering algorithms [111 
II 2\ . that are or will be included. All these 
new codes are written in the object-oriented 
programming language CH — h; in the beginning 
it was planned that at least the successors of 
PYTHIA and HERWIG would be based on the same 
underlying machinery and a repository of com- 
mon classes, CLHEP JT^j, which is also used in 
experiments. However, this turned out to be not 
feasible, and PYTHIA8 and SHERPA are relying on 
their own framework. 



2. PRESENTING SHERPA 

It is fair to state that, at the moment, SHERPA is 
the one of the new simulation tools, which is most 
advanced when it comes to the ability to actu- 
ally generating events. Currently, the following 
physics modules are implemented: 

• PDFs: 

Various PDFs - CTEQ [Hj and MRST Q2] 
in their original form as well as many other 
PDFs through the LHAPDF library QU - 
are interfaced. 

• Matrix elements: 

AMEGIC++ JJj is a matrix element gener- 
ator to describe hard scattering processes 
and decays at the tree level. Apart from the 
full SM AMEGIC++ contains the full MSSM 
[TB] in the notation of |19l2()j and an ADD 
model of large extra dimensions with 
its implementation described in [22j ■ SUSY 
particle spectra are provided through the 
SUSY Les Houches accord interface |23| . 

• Parton showers: 

For multiple QCD bremsstrahlung, i.e. 
the emission of secondary partons, 
SHERPA relies on APACIC++ which 
uses, similar to PYTHIA [7|, an ordering 
by virtuality supplemented with an explicit 
angular veto to ensure a proper treatment 
of quantum coherence. The merging of the 
hard matrix elements for multijet produc- 
tion and the subsequent parton shower is 
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achieved according to the merging proce- 
dure proposed in 25 26 and implemented 
in E7J. 

• Multiple parton interactions: 

A first simulation of the "hard" underly- 
ing event in the spirit of but supple- 
mented with parton showering has been im- 
plemented and tested. A more involved 
model is currently in preparation. 

• Hadronisation: 

The translation of the emerging partons 
into primordial hadrons is taken care of 
by the Lund string model |29|. This, as 
well as subsequent hadron decays are real- 
ized by an interface to the corresponding 
PYTHIA routines. However, a new version 
of cluster fragmentation j^O] is ready to be 
fully implemented in the near future. 

3. RESULTS 

3.1. Matrix elements 

AMEGIC++ automatically constructs Feynman di- 
agrams and helicity amplitudes [31132] for a given 
set of processes. For the helicity amplitudes, 
the formulation of |33| is employed. Having 
constructed them, AMEGIC++ simplifies and com- 
bines them by factoring out common parts and 
then writes the results out in library files to be 
compiled and linked with the core program. This 
leads to a drastic reduction of computing time. 
For the Monte-Carlo integration over phase space, 
the multi-channel approach of [34135] is being 
used. For each Feynman diagram, suitable phase 
space mappings are constructed and also writ- 
ten out as library files. During integration, the 
weight optimisation procedure selects successful 
channels having a large impact on the integra- 
tion. In AMEGIC++ , after a number of integration 
steps, these winning channels arc then further op- 
timised by employing VEGAS ^.3T| to select ran- 
dom numbers for them. 

AMEGIC++ has exhaustively been tested for a 
large number of production cross sections for six- 
body final states at an e + e~-collider [HZj and 
various processes at the LHC, see [35] ■ As 
an example for the latter, consider the pro- 



cesses pp — > e~v e + rcjets + X and pp — » 
e~v e bb + njets + X at the LHC. Cross sec- 
tions for these processes, obtained through 
Alpgen 1^1, CompHep gOI, MadEevent gT] 
14*2") . and AMEGIC++ can be found in Tablel3.ll 1 . 
In these comparisons, AMEGIC++ proved to work 
for up to eight external particles, and, thus, 
SHERPA includes a state-of-the-art matrix ele- 
ment generator, one of the key elements of mod- 
ern event generators. 

3.2. Merging of matrix elements and the 
parton shower 

In order to fully exploit the power of such multi- 
particle matrix elements, they have to be com- 
bined with the parton shower which models 
subsequent, secondary emission of softer QCD 
quanta. There arc different ways to do so, among 
them MC@NLO g3]. An alternative approach 
[25126) is to combine tree-level matrix elements 
for different jet multiplicities. This is done by 
defining two disjoint regions of jet production 
and evolution, separated by a jet measure de- 
fined according to the k± algorithm [44I45I46 1 ] . 
Then the matrix elements are reweighted with 
suitable Sudakov form factors such that the cor- 
responding matrix element becomes "exclusive" , 
and the parton showers are vetoed such that no 
extra jet is produced in the showering. This ap- 
proach guarantees independence of the jet separa- 
tion definition at leading logarithmic order. The 
algorithm has been implemented in full general- 
ity in SHERPA [22] 2 , forming one of its corner- 
stones. The approach and its implementation in 
SHERPA has been tested by comparing both with 
data and other codes, in a number of processes, 
cf. for example [51152153] , The findings in these 
comparisons were: 

• Self-consistency: 

Varying the k± cut in the internal jet defini- 
tion for the merging and the maximal num- 
ber of jets taken care off by the matrix el- 
ements, the approach has been found to be 
extremely stable and independent of the jet 

For details of the calculational setup and more results, 

cf. the MC4LHC homepage ESI- 

2 There exist some variations of this approach 47 48 49 50 
with different technical realisations of the same idea. 
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• Scale-independence: 

The shapes of characteristic distributions 
such as the transverse momentum of jets 
etc. are surprisingly stable (deviations of 
the order of 20% and less) under global vari- 
ations of the renormalisation and factorisa- 
tion scale in the matrix elements, Sudakov 
weights and the parton shower. The total 
cross section, being calculated at leading or- 
der only, however, depends much stronger 
on these choices. 



• Comparison with NLO results: 

The shapes of distributions obtained by 
SHERPA are in excellent agreement with 
those obtained by full NLO calculations. 

In this presentation, the merging approach will 
be applied to the case of jet production at Teva- 
tron, Runs I and II. At Run I, the DO collab- 
oration measured the ratio of the three-to-two 
jet rate R32, using the midpoint algorithm with 
different E± of the jets and in different regions 
of pseudorapidity of the jets [S^. The stability 
of the results obtained by SHERPA under varia- 
tions of renormalisation and factorisation scales 
Hr,f is demonstrated in Fig. In Fig. the 
SHERPA results are contrasted with those of a full 
NLO calculation In all these fig ures R32 is 
plotted against Ht, the scalar sum of the trans- 
verse momenta of all hard objects in the detec- 
tor. Jets are defined by E± > 40 GeV, \n\ < 3 
and R = 0.7. Again, stability of the results ob- 
tained by SHERPA under scale variations is found; 
the agreement with a full NLO calculation for the 
result is remarkable. These findings, however, are 
a continuation of what has been found before for 
other processes. 

At Run II, the DO collaboration measured the an- 
gular decorrelation in the azimuthal plane of the 
two leading jets in jet production |56| . In Fig. [3] 
the results of SHERPA are contrasted with the ex- 
perimental findings. Again, SHERPA is capable of 
precisely reproducing the QCD radiation pattern 
in these events. 
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Figure 1. R32 in dependence of Ht for different 
global factors on the renormalisation and factori- 
sation scale. 
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Figure 2. R32 in dependence of Ht as described 
by Sherpa and a full NLO calculation. 



4. CONCLUSIONS 

In this publication the event generator 
SHERPA has been presented. Results have been 
presented for the working of its internal matrix 
element generator AMEGIC++ , a state-of-the art 
tool. The implementation of the merging pro- 
cedure is a key ingredient of the SHERPA event 
generator. Results for jet production at the 
Tevatron prove that the merging of tree-level 
matrix elements and parton showers is work in 
a systematically correct manner. Further tests, 
especially in the simulation of more processes, 
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Figure 3. The azimuthal decorrelation of the two 
leading jets in the transversal plane. The four 
curves correspond to four different bins of p± of 
the leading jet. 



are ongoing. The results obtained so far indicate 
that SHERPA is perfectly suitable to meet the 
enhanced demands of the community to reliably 
simulate physics processes at the next generation 
of collider experiments. 
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